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Abstract 



We investigate one mechanism of the change in the isotopic composition of cosmologically distant clouds of 
interstellar gas whose matter was subjected only slightly to star formation processes. According to the standard 
cosmological model, the isotopic composition of the gas in such clouds was formed at the epoch of Big Bang 
nucleosynthesis and is determined only by the baryon density in the Universe. The dispersion in the available 
cloud composition observations exceeds the errors of individual measurements. This may indicate that there 
are mechanisms of the change in the composition of matter in the Universe after the completion of Big Bang 
nucleosynthesis. We have calculated the destruction and production rates of light isotopes (D, 3 He, 4 He) under 
the influence of photonuclear reactions triggered by the gamma-ray emission from active galactic nuclei (AGNs). 
We investigate the destruction and production of light elements depending on the spectral characteristics of 
the gamma-ray emission. We show that in comparison with previous works, taking into account the influence 
of spectral hardness on the photonuclear reaction rates can increase the characteristic radii of influence of the 
gamma-ray emission from AGNs by a factor of 2-8. The high gamma-ray luminosities of AGNs observed in recent 
years increase the previous estimates of the characteristic radii by two orders of magnitude. This may suggest 
that the influence of the emission from AGNs on the change in the composition of the medium in the immediate 
neighborhood (the host galaxy) has been underestimated. 
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1. Introduction 

The relative abundance of the light elements (H, D, 3 He, 
4 He, 6 Li, 7 Li) produced at the epoch of Big Bang nu- 
cleosynthesis is one of the criteria for testing the stan- 
dard cosmological model. According to this model, Big 
Bang nucleosynthesis occurred in the early Universe: at 
times from several minutes, when the Universe was ex- 
tremely hot (T> 10 9 K), to half an hour, when the tem- 
perature and density of the medium dropped to val- 
ues at which the nuclear reactions no longer proceeded. 
Subsequently, the composition of the interstellar and 
intergalactic medium changed mainly under the influ- 
ence of astration and radioactive decay. Therefore, by 
studying the isotopic composition of distant clouds of 
interstellar gas (z~2-=-3), where the star formation 
processes had not yet strongly affected the relative ele- 
mental abundances, their abundances at the completion 
time of Big Bang nucleosynthesis can be determined 

One of the most important aspects of such an anal- 
ysis is the abundance of deuterium. In comparison with 
other elements produced at the epoch of Big Bang nu- 
cleosynthesis, deuterium has the simplest nuclear struc- 
ture and is most sensitive to a change in the baryon 
density fif, - one of the key parameters of the standard 
cosmologicalmodel. The D/H ratio in the interstellar 
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medium of galaxies can be calculated from the absorp- 
tion by neutral gas of light emitted by distant quasars 
(QSO) by comparing the intensities of hydrogen absorp- 
tion lines with those of deuterium ones. When the pri- 
mordial D/H ratio is estimated, deuterium is commonly 
assumed to be only destroyed in the course of stellar evo- 
lution. Therefore, any D/H estimate is considered only 
as a lower limit for the primordial D/H ratio. 

Some of the first deuterium abundance measure- 
ments in clouds of interstellar gas at high redshifts gave 
large D/H ratios, ~ 1(T 4 (Webb et al. 1997; Rugers and 
Hogan 1996), which was inconsistent with the predic- 
tions of the standard model. An overview of the possi- 
ble mechanisms of the change in the abundances of light 
elements in clouds of interstellar gas relative to their 
primordial values was provided by Gnedin and Ostriker 
(1992). However, the new deuterium abundance obser- 
vations (Buries and Tytler 1998) disproving the previ- 
ous ones were in agreement with the theory, and the 
question about the change in abundance was essen- 
tially lifted. The currently available light-element abun- 
dance data taken from Steigman (2007) are presented 
in Table [TJ Eleven absorption systems of atomic deu- 
terium and hydrogen have been detected in the last ten 
years, and the D/H ratio has been measured in them 
with a high accuracy. The derived values are consis- 
tent, in order of magnitude, with the conclusions of 
the standard model, but there are some problems. At 
present, only seven systems in which the deuterium 
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abundance can be determined with good confidence 
from an analysis of absorption lines (Pettini et al. 2008, 
D/H=2.82±0.20 x 10~ 5 ) are used to estimate the D/H 
ratio. However, even for this sample the la scatter 
of D/H values was about 20%, which is considerably 
higher than the observational errors of individual mea- 
surements. 

It is believed that the scatter of deuterium abun- 
dances in such clouds can be provided by both star 
formation processes and deuterium depletion onto dust 
grains (Linsky et al. 2006). The interstellar gas that was 
subjected only slightly to chemical enrichment, as is sug- 
gested by the low relative abundance of heavy elements, 
O/H < 1/10 of the value measured in the Solar sys- 
tem, was studied in these seven systems. In the model 
of galactic chemical evolution, the decrease in the D/H 
ratio for such clouds is negligible, which was theoret- 
ically shown by Akerman et al. (2005), Prantzos and 
Ishimaru (2001), and Romano et al. (2006). Therefore, 
the large dispersion in the experimental data is most 
likely explained by an underestimation of the statisti- 
cal and systematic errors (Steigman 2007; Pettini et al. 
2008). Apart from the standard theory of Big Bang nu- 
cleosynthesis, the dispersion in the light-element abun- 
dance observations is associated with the evaporation 
of primordial black holes or the decay of long-lived X- 
particles. The high-energy hadrons and photons being 
produced during the decay of X-particles destroy the 
light elements synthesized in Big Bang nucleosynthesis. 
Effective He destruction and D, T, 3 He production can 
occur as a consequence of this process. These processes 
were investigated by Zeldovich et al. (1977), Kawasaki et 
al. (2005), and Carr et al. (2010). Stringent constraints 
were placed on the primordial abundance of the long- 
lived exotic X-particles. 

In this paper, we investigate one mechanism of the 
change in the relative composition of the medium that 
was previously considered by Boyd et al. (1989) and 
Gncdin and Ostriker (1992). Hard X-ray and gamma- 
ray emission is present in the spectrum of active galac- 
tic nuclei (AGNs). Propagating in the Universe, the 
high-energy gamma emission interact with the interstel- 
lar and intergalactic matter. The photodisintegration of 
4He with the production of D can be compared with the 
photodisintegration of D. This can cause both a decrease 
and an increase in the relative abundance of deuterium 
in some regions of interstellar gas. The magnitude of this 
effect depends on the initial composition of the medium 
(the abundance of 4 He in a cloud is greater than that of 
D by three orders of magnitude), the spectral hardness 
of the gamma-ray emission characterized by the param- 
eter r (see Eq. [3]), and its flux. 

Boyd et al. (1989) estimated this effect using the 
galactic center of NGC 4151 as an example. The au- 
thors showed that the gamma-ray source associated with 
the galactic center could cause 4 He photodisintegration 
and significant D production. However, the radius of 
influence of this process is very small and does not 
exceed 10 light-years in the case under consideration. 
The authors conclude that, given the rarity of gamma- 
ray sources in the Galaxy, the change in the isotopic 
composition of the Galactic medium is, in general, un- 
likely. Nevertheless, Casse et al. (1999) pointed out that 
blazars could affect the composition of the matter ac- 



Table 1. Primordial isotopic composition of matter in 
the Universe 



Element abundance 


Primordial value 6 


relative to hydrogen" 




= no/riH 


2.68xl0 -5 


-^3 He = n3 He/ n H 


1.06xl(T 5 


X 4 jj e = ni He /n H 


7.90xl(T 2 


X 7 n = m Li /n H 


4.30 xlO" 10 



a n - is the element number density in the medium. 
The data were taken from Steigman (2007). The abundances 
were reconciled with the baryon density Q.b obtained 
by analyzing the CMB anisotropy data. 



creting onto a black hole. When outflowing from the 
blazars central region, this matter can condense into in- 
tergalactic clouds. The authors also note that the final 
D abundance is sensitive to the spectral index of the 
quasar emission. However, the authors made no quanti- 
tative estimations. 

The goal of this paper is to consider the mechanism 
of element photodisintegration in a medium by the hard 
X-ray and gamma-ray emission from quasars by taking 
into account various spectral indices of emission and the 
high gamma-ray luminosities (compared to the values 
used by Boyd et al. 1989). 

2. THE CROSS SECTIONS FOR 
PHOTONUCLEAR REACTIONS 

One of the main points of our work is the determina- 
tion of photonuclear reaction rates. To calculate the re- 
action rates in a medium irradiated by a gamma-ray flux 
with various spectral indices T, the dependence of the 
reaction cross sections on the gamma-ray energy should 
be known. 

For deuterium (the simplest nuclear system), the de- 
pendence of the interaction cross section on the gamma- 
ray energy E can be calculated analytically (Bethe and 
Peierls 1935) and is written as: 



a(E) = a Q 



3/2 



(E-Q) 
E 3 



3/2 



(1) 



where ao is a dimensional constant, and Q is the reaction 
threshold. This dependence of the cross section on the 
gamma-ray energy describes well the experimental data 
on deuterium photodisintegration (see Fig. [T]). 

Other elements have a more complex nuclear struc- 
ture, and the theoretical calculations for many of them 
are not in satisfactory agreement with the experimental 
data (Dzhibuti et al. 1965; Balestraet al. 1977). Figure^ 
presents the experimental photonuclear reaction cross 
sections and their fits. The dependence of the reaction 
cross sections for other elements on the gamma-ray en- 
ergy is similar to the deuterium reaction cross section: 
there exist a threshold energy starting from which the 
reaction proceeds and a power-law decrease in the cross 
section. Taking this into account, we chose functions of 



Klimenko et al.: Influence of Gamma-Ray Emission 



3 




E [Me VI 



Figure 1. Cross section for the D(7, n)p reaction ver- 
sus 7-ray energy E. The dots indicate the experimental 
data from Hara et al. (2003), Birenbaum et al. (1985), 
Skopik et al. (1974), and Galey (1960). The solid curve 
was computed using Eq. ([1]) with 00 = 19.52mbarn, 
X 2 = 286.0. The dashed line represents fit ((2J with 
the following parameters: cto = 14.3 mbarn j3 = 1.19, 
7 = 2.55, x 2 = 173.5. 

the following form for the fit: 

a(E) = a ^- f L, (2) 

where ao, /3 and 7 are the parameters of the fit, and 
e = E/Q is the dimensionless energy. The best-fit pa- 
rameters are given in Table 2. For three 4 He destruc- 
tion reactions, the cross section near the threshold and 
at high energies cannot be described by one function. 
Therefore, for a more accurate fit, the energy range 
was divided into two parts, each of which was fitted 
independently. For the A He(j, 2p2n) reaction, there is 
a discrepancy in the absolute value of the cross sec- 
tion and the position of the maximum between the data 
from Arkatov et al. (1969, 1970), Gorbunov (1969), and 
Balestra et al. (1977). As was pointed out by Arkatov et 
al. (1970), this is because some arbitrariness is admitted 
when separating the individual events of the (^/,pn) and 
(7, 2p2n) reactions. In this case, the cross section was de- 
scribed by two functions, one of which was constructed 
from the data by Arkatov et al. (1970) and the other 
function was constructed from the data by Gorbunov 
(1969) and Balestra et al. (1977). However, the contri- 
bution from the 4 -ffe(7, 2p2n) reaction to the change in 
the abundance of the elements is so small that the differ- 
ence between the two fits does not lead to a noticeable 
change in the final result. 

Given the high gamma-ray energy, the photodisin- 
tegration reaction products can have a kinetic energy 
exceeding their binding energy. In this case, the thermal- 
ization of high-energy D and 3 He nuclei in the medium 
should be considered separately. Mostly the nuclei will 
collide with cold protons and will reduce its energy or 



be destroyed, depending on the energy of the colliding 
particles and the scattering pattern (elastic, inelastic). 
Deuterium has a low binding energy, Q^=2.23 MeV, 
will be most sensitive to the thermalization process. 
The width of the reaction cross section maximum is 
A E ~ 2 10 MeV (see Fig. CU and Table The ki- 
netic energy of the produced nuclei can be estimated as 
Ek — A#/2, which can exceed the deuterium binding 
energy Qd. Therefore, the thermalization of D and 3 He 
nuclei in themedium can partially destroy the produced 
amount of elements. 

Thus, complex statistical calculations of the ther- 
malization of nuclei required to reliably estimate the 
amount of produced D and 3 He nuclei. In this paper, 
we consider the effect that consists in a great change in 
the abundance of light nuclei in the medium irradiated 
by emission with various values of the spectral hardness 
r. This effect turns out to be especially important in 
the immediate vicinity of AGNs. 



3. OBSERVATIONAL DATA ON ACTIVE 
GALACTIC NUCLEI 

Blazars and quasars (AGNs) are extremely intense 
sources of energy release in the Universe; they also have 
enormous gamma-ray luminosities. Since gamma rays 
with an energy above 2.23 MeV (the deuterium binding 
energy) are involved in the photodisintegration of ele- 
ments, we are interested in the gamma-ray energy range 
from 1 to 1000 MeV. The upper boundary was chosen 
with some margin: since the quasar spectrum falls off as 
a power law, taking into account the energy range from 
100 to 1000 MeV causes the reaction rates to increase 
by no more than 1%. 

The spectra of such objects were investigated at the 
orbital observatories: AGILE (30-500 MeV; see, e.g., 
D'Ammando et al. 2011), COMPTEL (0.75-30 MeV; 
see, e.g., Schonfelder et al. 2000), EGRET (30 MeV- 
20 GeV; see, e.g., Hartman et al. 1999), and Fermi LAT 
(20 MeV-300 GeV; see, e.g., Abdo et al. 2010). 

In the region of photonuclear reaction thresholds (2- 
30 MeV) , there are data on the spectra from the survey 
by Schonfelder et al. (2000). As was shown in the surveys 
by Hartman et al. (1999) and Schonfelder et al. (2000), 
the quasar emission at high energies (E 1 > 1 MeV) is 
well fitted by a power law: 

F^(E)=Fo(^) 1 (3) 

where Fo is a flux density in phot cm 2 s -1 MeV -1 , T 
is the spectral index, and £0 ~~ is the scale factor. Their 
typical values are F ~ 10 -8 phot cm 2 s _1 MeV -1 , T = 
1 4, £ = 1 + 300 MeV. 



4. KINETICS OF THE ISOTOPIC 

COMPOSITION OF THE INTERSTELLAR 
MEDIUM 

The kinetics of the isotopic composition of the inter- 
stellar medium as a function of the exposure time and 
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Figure 2. Photonuclear reaction cross sections versus 7-ray energy E. The solid curves represent fit ^ with 
<7o,/3,7 from Table [5J The dots indicate the experimental data: T( 7 , n)D Varfolomeev and Gorbunov (1965), 
Bosch et al. (1964), Skopik et al. (1981), Faul et al. (1981); T( 7 ,2n)p Faul et al. (1980). 
3 He(7,p)D Ticcioni et al. (1973), Stewart et al. (1965), Warren et al. (1963), Naito et al. (2006); 
3 He(7,n)2p Fetisov et al. (1965), Ticcioni et al. (1973), Faul et al. (1981), Naito et al. (2006); 
4 He(7,p)T Shima et al. (2001,2005), Kusakabe et al. (2009), Arkatov et al. (1978b), Gorbunov (1968), 
Balestra et al. (1977); 4 He( 7 ,n) 3 He Irish et al. (1975), Balestra et al. (1977), Gorbunov (1968), 
Shima et al. (2001,2005), Kusakabe et al. (2009), Nilson (2005,2007), Arkatov (1978a); 
4 He( 7 , D)D Arkatov et al. (1980); 4 He(7,pn)D Shima et al. (2001,2005), Arkatov et al. (1969b), Gorbunov (1958); 
4 He(7,2p2n) Balestra et al. (1977), Arkatov et al. (1969a), Gorbunov (1969). 



spectral parameters can be obtained by solving the sys- 
tem of ordinary differential equations (|U - [7]) : 



He 



dt 



4 He 



dX T 



dt 



~ {kr^He 

I- ki He ^3 H Xi He 
dXs He 



t-p) Xi 



dt 



= -{k 



3 He^D 



^He 



_2p) X3 H 



+ k4 He _ > 3 He X4 He + k T ^3 He X T 

dX D 



dt 



-kD^p+nXu + (2 • kiHe^D+D 



k±He^D+pn)XiHe + k3 He ^ D Xs He + kx^D^T, 



(4) 
(5) 
(6) 
(7) 



where X denote the elemental abundances, and k are the 
photonuclear reaction rates. The element abundances 
calculated using the theory of Big Bang nucleosynthesis 



with the baryon density £Ib obtained by analyzing the 
cosmic microwave background (CMB) anisotropy data 
(see Table [I| were taken as the medium's initial com- 
position. The influence of the elements with an atomic 
weight greater than that of 4 He on the change in the 
abundances of lighter elements may be neglected, be- 
cause their relative abundance in a medium similar to 
the primordial one is low (( 7 Li/H)< 5 x 10 -10 ), while 
the cross sections are of the same order of magnitude 
(see Berman et al. 1965; Varlamov et al. 1986). Using 
Eqs. (J2J) and ([3]), the photonuclear reaction rates can be 
written as: 



k = J a (E) F 7 (E) dE = 10~ 27 (- r " 



x F e 



Q 



1 MeV 



1-r 



1000/Q 



mbarn 





(8) 



,7+r 



■ dx, 
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Table 2. Parameters of the fit to the photonuclear cross sections . 
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Reaction 


Q, MeV 


E7, MeV 


o"o, mbarn 


P 


7 


D(7,n)p 


2.23 




14.3 


1.19 


2.55 


T(7,n)D 


6.257 




9.86 


1.7 


3.6 


T(7,2n)p 


8.482 




22.1 


2.11 


4.65 


He{-y,p)U 






Ol A 


O A Q 

z.4o 


A 


3 He(7, n)2p 


7.72 




9.48 


1.87 


3.43 


4 Hef7,p)T 


19.81 




E-, < 40 
E-, > 40 


145.5 
19.5 


2.3 
1.0 


7.4 
4.5 


4 He(7,n) 3 He 


20.57 




E-y < 50 
F 7 > 50 


60.2 
9.11 


1.8 
0.6 


6.1 
3.5 


4 He( 7 ,D)D 


25.0 




F 7 < 50 
F 7 > 50 


0.64 
0.01 


1.95 
0.1 


10.0 
4.3 


4 He(7,pn)D 


28.0 




2.1 


1.5 


3.6 


4 He( 7 , 2p2n) 


28.0 




0.81 


1.1 


3.5 



a The fitting formulas give an accuracy of at least 10%. 

Note. The spontaneous decay rate of tritium T(/3") 3 He is k = 1.82709 ± 0.0026 x 10" 9 s _1 , see Akulov and Mamyrin (2004). 
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Figure 3. Effective D, 3 He and 4 He production and destruction rates when primordial matter is exposed to 
7-rays with energy E = 1-1000 MeV with the spectrum dN/dE = Nq(Eq/E) t versus T. Nq corresponds to the 
7-ray flux with L 7 = 10 47 erg s^ 1 at the distance of 1 kpc. 



where the dependence on /3, 7 and T is determined by 
an integral that can be expressed in terms of Eulers 
incomplete beta functions: 



1000/Q 



(x-lf 

y 7+ r dx = Bi(r + 7 - - 1,0 + 1) + 



+ B Q/100 o(r + 7-/3-l,/3 + l). 



(9) 



The system of equations (JU)-© is solved analytically. 
In view of the fitting formulas © and the element 
abundances at any instant of time are expressed in alge- 
braic form in terms of special functions. This allows the 
distributions of elemental abundances in the medium 
to be rapidly calculated for various spectral parame- 



ters. The results of our calculations are presented in 
Figures SHHl 

The change of the element abundances for a fixed 
spectral slope T is determined by only one parameter 
x = Foept, the multiplication of the gamma-ray flux 
at an energy of 1 MeV by the exposure time. We as- 
sume that the gamma-ray flux is produced by the cen- 
tral source and is determined only by the distance to 
the quasar and its luminosity. Taking into account the 
correlation between the optical and gamma-ray lumi- 
nosities (Arshakian 2011), we may consider the opti- 
cal Eddington luminosity \jEdd — ^nGMm p c/ax — 
1.5 x 10 46 Mg erg s _1 , where ctt is the Thomson cross 
section, as an upper limit for the gamma-ray luminosity 
of sources. The range of admissible cloud exposure times 
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Figure 4. Total luminosity of the 7-ray source L 7 in the energy range 1-10 MeV versus spectral slope V 
corresponding to a ±10% change in the abundance of D and 3 He in the primordial matter at a distance of 1 kpc 
in an exposure time of 10 9 yr. 



2,0 



1.5 



1,0 - 



0,5 



0,0 





■] n 1 


1 ji 


I 1 




. — 1.5 






2.0 


4 

1 




- — 3.1c 






— 3.6; 


1 




I 


; 




■ a / 


1 






f 

' a 


r 


v*o.y 


" H * 


1 I 

^ 1 




j? 


/ 1 






_ ^ I 




V-Q.1 






1 ■! 1 







10' 




F Q s [photons/ cm s MeV] F Q s [photons/ cm s MeV] 



Figure 5. Change of the D (a) and 3 He (b) abundances in a medium exposed for 10 9 yr versus flux density at 
an energy of 1 MeV FqEq. As the spectral slope T increases, the curve is shifted rightward; at V > = 3.13, 
the Ad curve acquires a negative value (D is destroyed). T — 3.62 corresponds to the minimum flux density 



FqEq 



10 11 phot cm 2 s 1 MeV 1 at which 10% of D will be destroyed in 10 9 yr. 



is estimated via the lifetime of a quasar at its active 
phase t= 10 5 -f- 10 9 yr (Martini and Weinberg 2001). We 
consider the characteristic exposure time to = 10 9 yr to 
be able to compare our results with the work by Boyd et 
al. (1989) and to take into account the maximally pos- 
sible effects when estimating the deuterium production. 

For the limiting gamma-ray luminosity L 7 = 
10 47 erg s _1 (i.e., the Eddington luminosity of a quasar 
with a mass of 10 9 M sun ), the element abundances 
changes by ~ 10% at distances of ~ 100 pc. The self- 
similarity of the model allows the radius of influence of 
the quasar (on the medium's composition) to be recal- 
culated for such quasars of a different luminosity and 
for a different exposure time by rescaling: 

/ T f \ (1/2) 

' (10) 

In contrast to previous works, we also investigated 
the dependence of the reaction rates on the spectral 



index. For a quasar with a fixed luminosity L 7 = 
10 47 erg s _1 in the energy range 1-1000 MeV, the D 
and 3 He production and destruction rates are plotted 
against the spectral hardness T in Fig. [3J These plots 
describe the dynamics of the medium's composition on 
short exposure time scales, when the relationship be- 
tween the abundances of 4 He and 3 He, D in the medium 
did not change significantly. The maximum in the reac- 
tion rates (at r ~ 2) is determined by the relationship 
between the position of the maximum in the reaction 
cross sections and the power-law shape of the spectrum. 
For deuterium (Fig. [5b), there is a feature T — 3.1. For 
harder spectra (r < 3.1), deuterium in the medium will 
be only produced (through the photodisintegration of 
4 He); conversely, for softer spectra (r > 3.1), the de- 
struction of deuterium can exceed its production. There 
are no such features for 3 He (see Fig. [3^). it is only pro- 
duced. The contributions from the 3 He production dur- 
ing the 4 He destruction exceeds the 3 He photodisinte- 
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gration by more than two orders of magnitude. Since 
there are no 4 He sources in the model under consid- 
eration, 4 He is only destroyed. The amount of tritium 
is established at a constant level (until the 4 He abun- 
dance changes significantly) determined by the equilib- 
rium condition dXT/dt = 0: 



X 



ki 



4 He- 



He^T 



K 7 



Xi He ({)) kiH ^ T exp(-if4 ge f), 



(11) 

where Ki He and Kt denote the total element destruc- 
tion rates. 

The features in the relationships between the reac- 
tion rates determine the pattern of changes in the ele- 
ment abundances. Figure 0] define the regions of spectral 
parameters (the gamma-ray luminosity of the quasar L 7 
for an isotropic source and the spectral index V) corre- 
sponding to a 10% change of the elemental abundances 
in a cloud Apj and A 3 jj e at a distance of 1 kpc in an 
exposure time of 1 Gyr. As has been pointed out above, 
the destruction of D is possible in a certain domain of 
parameters. In Fig. [SI the change in the abundances of D 
and 3 He is plotted against the spectral flux density F Eq 
at an energy of 1 MeV for various spectral indices T. The 
abundance of D and 3 He is determined only by the ex- 
posure parameter x oc Fo ■ t. Therefore, apart from the 
spatial change in elemental abundances (the flux den- 
sity decreases as Fo oc r~ 2 ), these plots illustrate the 
temporal evolution of the D and 3 He abundances in the 
medium. 

Another feature in the distribution of deuterium is 
that when the medium is irradiated by emission with a 
soft spectrum, the decrease in the mass fraction of D is 
subsequently (at high flux densities or on long exposure 
time scales) replaced by a significant increase. This is 
because a new D production channel appears from the 
newly produced 3 He. 



5. THE EFFECTIVE RADIUS OF INFLUENCE. 
AN ISOTROPIC SOURCE 

The effective radius of influence of a quasar as a function 
of its gamma-ray luminosity and spectral hardness in a 
time of ~1 Gyr can be determined within the spherically 
symmetric quasar emission model. As has been shown 
above, the characteristic radius of influence of a quasar 
with a luminosity L 7 = 10 47 erg s _1 in the energy range 
1-1000 MeV is ~ 100 pc. The Hubble expansion may be 
neglected at such distances, and the flux density is then 
related to the luminosity by the standard relation 



L-, 



4irR 2 



1000 

J F^(E)EdE = F ell(T), (12) 



the latter expression was derived by taking into account 
Eq. ©, where I(T) = 100 f'^ 1 and H T ) = m(1000) at 

r = 2.0. 

Figures [5] and [7] show the profiles of the D/H and 
3 Hc/H ratios as a function of the quasar distance. The 
distribution of elements in the medium changes greatly 
as the quasar is approached. The main change occurs at 
^100 pc, where the amount of D and 3 He increases by 
two or three orders of magnitude. At shorter distances, 



effective photodisintegration of both 4 He and D and 3 Hc 
takes place. Since there are no He sources in the model 
under consideration, the He abundance decreases ex- 
ponentially with increasing flux density or decreasing 
quasar distance. Thus, the effect of a change in the D 
and 3 He abundances is strongly local (within ~1 kpc of 
the AGN center). 

The changes in the element abundances are also sen- 
sitive to the spectral hardness T. For different spectra, 
the elemental abundances in the medium can differ by 
an order of magnitude at the same exposure time and 
quasar luminosity. Figures [6] and [7] present several curves 
corresponding to different spectral slopes. The depen- 
dence of the radius of influence on the spectral index T 
is determined by the dependence of the reaction rates 
(see Fig. [3]): D and 3 He are produced in large quantities 
and at great distances at T ~ 1.8 and 2.0, respectively. 



6. MOTION OF THE MEDIUM 

The distributions of the abundances of light elements 
formed through photonuclear reactions (see Figs. H2 and 
[7]) were obtained by assuming the medium to be at rest, 
when the matter at different distances from the quasar 
center is not mixed. In reality, the matter in the in- 
terstellar medium moves, and the distributions of the 
elements change. 

Estimating the global effect of the contribution from 
the emission of quasars to the change of the D and 3 He 
abundances in the Universe as a whole (after the com- 
pletion of Big Bang nucleosynthesis), we perform a uni- 
form averaging of the changed D and 3 He abundances 
over the space in which the change occurred: 



< X > (R,L 7 ,T) = 



n _ 1 J R X(r)4r 



V 



(13) 



where X denotes the abundances of the elements: (D/H), 
( 3 He/H) and ( 4 He/H). 

The upper limit of the integral R in Eq. (|T5|) . which 
we call the effective zone of influence of a quasar, is 
chosen from the conditions under which the quasar in- 
fluence will change the composition of the medium by 
10% at 1 Gyr. The limiting quasar luminosity L 7 
1 1 1 erg s _1 and exposure time t=10 9 yr are used to es- 
timate the maximum influence of the quasar emission on 
the composition of the medium. The size of the region 
of influence is ~1 kpc for D and ~ 10 kpc for 3 He. The 
deviation of the averaged abundance of D and 3 He from 
its primordial value does not exceed 5%; the produced D 
and 3 He are plotted against the quasar distance for var- 
ious spectral slopes in Figs. H] and O The element abun- 
dance averaged over a larger volume decreases in inverse 
proportion to the distance squared, < X >oc R~ 2 . 

Thus, the global effect of the gamma-ray emis- 
sion from quasars (at a limiting isotropic luminosity of 
10 47 erg s _1 ) on the composition of the medium turns 
out to be insignificant, and the emission from quasars 
cannot be considered to be a source of the observed 
excesses of the D and 3 He mass fractions above the pri- 
mordial values. 

However, this estimate obtained by assuming uni- 
form mixing is only a lower limit for the influence of a 
quasar on the change in the isotopic composition of the 
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R, [kpc] 

Figure 6. D/H ratios in the medium versus distance to a quasar of luminosity L 7 = 10 47 erg s _1 . The 
D abundance relative to its primordial value (D/H) pr = 2.68 x 10~ 5 is along the vertical axis. The curves are 
presented for various spectral slopes T = 1.0, 1.5, 1.8, 2.7, 3.8 of a source with the same luminosity. V = 2.7 
corresponds to the spectral slope for the quasar emission considered by Boyd et al. (1989). The largest amount of 
D is produced at T = 1.83. 
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Figure 7. 3 He/H ratio in the medium versus distance to a quasar of luminosity L 7 = 10 47 erg s _1 . The 3 He 
abundance relative to its primordial value ( 3 He/H) pr = 1.06 x 10 -5 is along the vertical axis. The curves are 
presented for various spectral slopes T= 1.0, 2.0, 2.7, 3.6 of a source with the same luminosity. T = 2.7 corresponds 
to the spectral slope of the quasar emission considered by Boyd et al. (1989). T = 2.0, the maximum amount of 
3 He is produced in the medium. 




surrounding medium. Considering the model of nonuni- 
form mixing, a situation is possible where a cloud of gas 
overrich in D and 3 He can be accelerated by a jet and 
be ejected into the medium surrounding the quasar. The 
matter from the interstellar medium can also fall to the 
quasar through an accretion disk, be exposed, and then 
ejected back. The existence of such a mechanism was 
previously pointed out by Casse et al. (1999). 

The pollution of the medium through this mecha- 
nism depends on many parameters, such as the relation- 
ship between the required matter exposure time near the 



quasar and the cloud lifetime near the galactic nucleus, 
the degree of destruction of the produced elements dur- 
ing the ejection of matter into the intergalactic medium, 
and the condensation of the ejected matter into small in- 
tergalactic clouds. This mechanism seems possible, but 
a detailed study of the theory of accretion disks around 
AGNs, the theory of jets, and other theories is needed 
for its assessment. 
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Figure 8. Average D/H ratio versus size of the averaging region corresponding to a change in the composition 
of primordial matter when exposed to 7-rays with an energy of 11000 MeV from a source of luminosity L 7 = 
10 47 erg s _1 . D/H in units of its primordial value (D/H) pr is along the vertical axis. The curves are presented for 
several spectral slopes r=1.0, 1.5, 1.83, 2.7, 3.8. T — 1.83 corresponds to the maximum amount of D produced in 
the medium. Boyd et al. (1989) used the model of a source with a spectral slope T = 2.7. 




R, [kpc] 

Figure 9. Average ( 3 He/H) ratio versus size of the averaging region corresponding to a change in the composition 
of primordial matter when exposed to 7-rays with an energy of 11000 MeV from a source of luminosity L 7 = 
10 47 erg s _1 . 3 He/H in units of its primordial value ( 3 He/H) pr is along the vertical axis. The curves are presented 
for several spectral slopes V = 1.0, 1.5, 2.0, 2.7, 3.6. T = 2.0 corresponds to the maximum amount of 3 He produced 
in the medium. Boyd et al. (1989) used the model of a source with a spectral slope T — 2.7. 



7. DISCUSSION 

Previously, the influence of photodisintegration pro- 
cesses on the abundances of light elements was con- 
sidered by Boyd et al. (1989) using the galactic cen- 
ter of NGC 4151 as an example for the chosen spectral 
slope r = 2.7. The gamma-ray flux with E > 2 MeV 
was assumed to be F = 4 x 10 16 phot cm 2 s _1 at 



a distance of 2 lt-days (the corresponding luminosity 
is L 7 ~ 1.6 x 10 44 erg s _1 in the energy range 1— 
1000 MeV). On the whole, our calculations reproduce 
the results of this paper, but there are also differences 
(see Fig. [TO]) . In particular, to estimate the reaction 
rates, the authors used reaction cross sections (nor- 
malized to the gamma-ray flux) that turned out to be 
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slightly overestimated. Table [3] compares the data from 
Boyd et al. (1989) with the values obtained by our fit- 
ting of the new data. The difference between the cross 
sections for some reactions reaches 200%. 

Using the same reaction cross sections and initial ele- 
ment mass fractions as those taken by Boyd et al. (1989), 
we reproduced the results of this paper. However, in 
the region r<10 lt-years, where the authors disregarded 
the change in the abundance of 4 He, the results dis- 
agree. According to our calculations, as the galactic cen- 
ter is approached, the amount of D and 3 He does not 
reach a constant level but reaches its maximum value 
at r ~ 4 — 7 light years. Subsequently, as the galac- 
tic center is approached, the element abundances will 
not remain constant, because the photodisintcgration 
of all elements begins to dominate as the flux increases. 
Accordingly, the change in the abundance of 4 He should 
be taken into account to correctly estimate the influ- 
ence of gamma-ray emission on the isotopic composi- 
tion of the medium in the immediate neighborhood of 
the AGN. 

The depletion of 4 He was taken into account by 
Balbcs et al. (1996), who also considered the influence 
of photonuclear reactions on the primordial composition 
of the medium in the Universe. Figure [T2l compares the 
results of our calculations with those from Balbes et al. 
(1996). Our calculations reproduce the results of this 
paper. Nevertheless, these authors did not consider the 
dependence of the D and 3 Hc production on the spec- 
tral parameters of the emission L 1 and T by assuming 
that the cross sections for spectra of different hardness 
increased by no more than a factor of 4, which actually 
is not the case. When varying the spectral hardness, the 
reaction cross sections change by one or two orders of 
magnitude. 

In contrast to Boyd et al. (1989) and 
Balbes et al. (1996), we considered the influence 
of spectral hardness T on the produced amount of D 
and 3 He. We found that the amount of D and 3 He 
produced in the zone of influence of a quasar for spectra 
of different values of the spectral index could exceed 
the value found by Boyd et al. (1989) by a factor of 7 
for D and a factor of 3 for 3 He. We retained the quasar 
luminosity L 1 — 1.6 x 10 44 erg s _1 , changing only the 
distribution of gamma-ray photons over the spectrum 
(various spectral indices). The derived dependencies 
are shown in Fig. [TT] Thus, the effect studied by Boyd 
et al. (1989) using the galactic center of NGC 4151 as 
an example may turn out to be negligible if soft spectra 
with r > 3 are considered or be more significant for 
harder spectra with r ~ 2 relative to T — 2.7 used by 
Boyd et al. (1989). 

New observations of quasars and blazars (see, e.g., 
Arshakian 2011) give an upper limit on the blazar lumi- 
nosity in the energy range 1-1000 MeV of 10 49 erg s _1 . 
This increases the effective zone of influence by more 
than three orders of magnitude (provided that the ac- 
tivity time of the nucleus is the same). Note, however, 
that the emission for blazars is localized in a narrow 
region corresponding to the jet and only in this region 
is a high value of the characteristic distance of the in- 
fluence of emission on the medium possible. According 
to our calculations, it is ~1 kpc for D and ~10 kpc 
for 3 He. In the case where the isotropic source emission 



Table 3. Comparison of the reaction cross sections. 



< a >, mbarn 



Reaction 


Boyd et al. (1989) 


our calculation 


D(7, n)p 


1.62 


1.33 


3 He( 7 ,p)D 


0.103 


0.073 


3 He(7,n)2p 


0.066 


0.050 


4 He( 7 ,p)T 


0.017 


0.015 


4 He(7,n) 3 He 


0.017 


0.015 


4 He( 7 ,D)D 


5.5xl0 -5 


2.62xl0~ 5 


4 He(7,pn)D 


0.0013 


0.00124 



model can be applied, the gamma-ray luminosity of a 
quasar typically does not exceed 10 47 erg s _1 , i.e., the 
optical Eddington luminosity of a quasar with a mass 
of 10 9 M sun . In this case, the distance at which the 
D and 3 He abundances change locally are, respectively, 
^100 pc and 1 kpc for 3 He. The derived radii of influ- 
ence turned out to be comparable to the characteristic 
distances at which the structure of the quasar itself (the 
source being non-point like, the accretion disk, and the 
jet) should be taken into account. A detailed account of 
the objects structure can cause the effects to increase. 
In the process of its accretion onto an active nucleus 
through an accretion disk, being in the immediate vicin- 
ity of the active nucleus, the matter can be irradiated by 
considerably larger gamma-ray fluxes and, accordingly, 
can change its isotopic composition in a shorter time. 
The D- and 3 He-overrich matter ejected into the inter- 
stellar and intergalactic medium can cause the isotopic 
composition of the medium at great distances from the 
AGN to change. 

As an example, the application of the examined 
mechanism of the change in the isotopic composition 
of the medium to our Galaxy can be considered for 
the so-called ' 3 He Problem'. Thus, the absence of a 
3 Hc/H gradient (Bania et al. 2007) in our Galaxy al- 
lows a limit to be placed on the activity time of the 
Galactic nucleus in the past. The hyperfine transition 
line of 3 He+ (8.665 GHz) observed in HII regions and 
planetary nebulae is used to determine the amount of 
3 He in the interstellar medium. The 3 He/H ratio is de- 
termined via the 3 He+/H+ ratio. The 3 He abundance 
at Galactocentric distances from to 16 kpc turns out 
to be at the same level, 3 He/H=1.79±0.65 x 10" 5 . At a 
smaller Galactocentric distance, 3 He was observed only 
in the HII region of the source Sgr B2 in 1990 (Balser 
1994): He 3 /H=2.49xl0~ 5 , the distance to Sgr B2 was 
found to be r=0.09±0.03 kpc (Reid et al. 2009). 

Assuming that the increase in 3 He/H relative to its 
primordial value ( 3 He/H) pr was caused only by pho- 
tonuclear reactions, we can estimate the activity time 
of the Galactic nucleus in the past. The mass of the 
supermassive black hole at the center of our Galaxy, 
M BH = 4.5 ± 0.4 x 10 6 M sun (Ghez et al. 2008), turns 
out to be so small that to produce the observed amount 
of 3 He at a distance of 100 pc, the medium should 
be exposed at the Eddington luminosity limit of the 
Galactic nucleus l^Edd — 6.75 x 10 44 erg s _1 for to < 
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Figure 10. ( 4 He/H), ( 3 He/H) and (D/H) abundances versus distance to the galactic center of NGC 4151. The 
circles indicate the 3 He/H and D/H curves from Boyd et al. (1989). The solid curves represent our calculations 
with the data from Boyd et al. (1989). The dashed curves represent our complete calculations. 
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Figure 11. Relative amount of produced D and 3 He through photonuclear reactions in the zone of influence of 
the galactic center of NGC 4151 for 1 Gyr versus spectral hardness T. The values are normalized to the amount 
of D and 3 He produced when the medium is exposed to emission with a spectrum with T = 2.7 (this value was 
considered by Boyd et al. (1989)). The maximum amount of D is produced at T = 1.83 and exceeds the value from 
Boyd et al. (1989) by more than a factor of 7. The maximum amount of 3 He corresponds to T = 2.0. 



1.7 ± 0.9 x 10 9 yr (see Fig. IT3"|) . which exceeds consid- 
erably any reasonable estimates of the lifetime for the 
nuclei. Some indirect data on the deuterium abundance 
point to the existence of a gradient with Galactocentric 
distance (Lubowich and Pasachoff 2009). This can be 
due to both deuterium astration in stars and complex 
chemistry of the interstellar medium. 

8. CONCLUSIONS 

We considered the influence of photonuclear reactions 
on the change in the composition of the interstellar 



medium. Directional exposure of the medium to gamma- 
ray emission was shown to cause the relative D and 3 Hc 
abundances in the medium to increase. 

We calculated the photonuclear reaction rates as a 
function of the spectral hardness of the emission. We 
fitted the photodisintcgration reaction cross sections by 
analytical formulas using which the reaction rates can 
also be expressed analytically in terms of Eulers beta 
functions. 

In comparison with previous works (Boyd et al. 1989; 
Balbes et al. 1996; and others, L = 10 44 erg s _1 , T = 2.7) 
, we considered power-law spectra with various values of 
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Figure 12. ( 4 He/H), ( 3 He/H) and (D/H) abundances versus exposure parameter x = Ljt/Awr 2 . The dots indicate 
the 4 He/H, 3 He/H, and D/H curves from Balbes et al. (1996) for ijio = 100 (Fig. 1). The solid curves represent 
our calculations. 
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Figure 13. The dots indicate the experimental data on the 3 He abundance in the Galaxy as a function of 
the Galactocentric distance. The nearest measurement was made in the HII region of Sgr B2. The dashed curve 
indicates the theoretical calculation of the 3 He abundance in a region exposed to gamma-ray emission from the 
galactic center at the Eddington luminosity limit for 10 9 yr. The 3He abundance curve for spectra of different 
hardness is within the measurement accuracy of the distance to Sgr B2. 



r = l-j-4 and took into account the high gamma-ray lu- 
minosities of quasars, L~ 10 47 erg s^ 1 . Allowance for 
the spectral hardness in the model under consideration 
leads to a change in the final mass fractions of the ele- 
ments by a factor of 37. 

The global effect from the production of D and 3 He 
through the quasar emission in the Universe, on the 
whole, turns out to be insignificant. However, this is not 
the case if the changes in the composition of the medium 
are considered in a local region. The ejected matter with 
a changed isotopic composition can condense into clouds 



of interstellar and intergalactic gas, thereby changing 
the composition of such clouds. A further development 
of the model requires a detailed analysis of the structure 
of AGNs, i.e., allowance for the extent of the jets, the 
finite sizes of the accretion disk, and the possibility of 
the matter processed in the disk being ejected into the 
interstellar medium, which can cause the isotopic com- 
position of the medium at great distances to change. 

Within the model under consideration, the absence 
of a 3He/H gradient in the Galaxy constrains the activ- 
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ity time of the Galactic center at its Eddington luminos- 
ity limit (L= e^xlO^ergs- 1 )^ < 1.7 ±0.9 x 10 9 yr. 
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